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Abstract 

We conduct a study based on the Bures (minimal monotone) metric, analogous to that recently 
reported for the Hilbert-Schmidt (flat or Euclidean) metric (arXiv:0704.3723v2). Among the 
interesting results obtained there had been proportionalities — in exact correspondence to the 
Dyson indices (3 = 1,2,4 of random matrix theory — between the fourth, second and first powers 
of the separability functions St ype {p) for real, complex and quaternionic qubit-qubit scenarios, 



Here u = y pUpH » with p being a 4 x 4 density matrix. Separability functions have proved 
useful — in the framework of the Bloore (correlation coefficient/off-diagonal scaling) parame- 
terization of density matrices — for the calculation of separability probabilities. We find — for 
certain, basic simple scenarios (in which the diagonal entries of p are unrestricted, and one 
or two off-diagonal [real, complex or quaternionic] pairs of entries are nonzero) — that these 
proportionalities no longer strictly hold in the Bures case, but do come remarkably close to holding. 
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In our recent study [1], we reported a number of developments of value in resolving the 
clearly challenging and conceptually important question of the probabilities — using the mea- 
sure induced by the Hilbert-Schmdit (HS) metric — that generic real, complex and quater- 
nionic qubit-qubit and qubit-qutrit states are separable/disentangled. An essential com- 
ponent in this progress was the use of the (simple) parameterization of the n x n density 
matrices (p) that had been originally proposed by Bloore Q . This involves reparameterizing 
the off-diagonal entries pij as y/paPjj (xij + Wij)- In the real case (y^ = 0), the nxn matrix of 



Xij% being necessarily nonnegative definite, with Xij G [- 



ll, has the form of a correlation 



|5|]. (Correlation matrices can 



n 

matrix — a basic object of study in descriptive statistics [3, 
be obtained by standardizing covariance matrices. Density matrices have been viewed as 
covariance matrices of multivariate normal [Gaussian] distributions [6] . The possible states 
of polarization of a two-photon system are describable by six Stokes parameters and a 3 x 3 
"polarization correlation" matrix [?]].) 

A major virtue (of course, unrecognized more than thirty years ago in the 1976 paper 
) of the Bloore parameterization is that it allows one to express the well-known Peres- 



Horodecki positive-partial-transposition criterion [8|, [9], [h]] for separability using fewer vari- 
ables than one would naively anticipate. Since we are largely concerned with the evaluation 
of high- dimensional integrals, this reduction in number of relevant variables is certainly of 
considerable importance. 

Here, we parallel the sequential approach of Zyczkowski and Sommers in that they, first, 
computed the total volume of (separable and nonseparable] nxn density matrices in terms 
of the (flat or Euclidean) Hilbert-Schmidt metric llj 12j, sees. 9.6-9.6, 14.3], and then, 
using the fundamentally important Bures (minimal monotone) metric 12, sec. 14.4] l^ . 
(In particular, they employed the Laguerre ensemble of random matrix theory 14J in both 
sets of computations (cf. 15]). The Bures and HS metrics were compared by Hall l^J, who 
concluded that the Bures induced the "minimal-knowledge ensemble" (cf. [17]).) That is, 
we will seek now to extend the form of analysis applied in the Hilbert-Schmidt context in 
tlfl to the Bures setting. 

n 

To begin, let us review the most elementary findings reported in [l|, sec. II.A.l]. The 
simplest (four-parameter) scenario studied there posits a 4 x 4 density matrix p with fully 
general diagonal entries (pu, P22, P33, Pu = 1 — Pu — P22 — P33) and only one pair of real 
off-diagonal non-zero entries, P23 = P32- The HS separability function for that scenario was 



2 



found to take the form 



t, eq. 



(20)], 

,2a < a < 1 

> 1 

where we will now primarily employ (purely as a matter of convenience) the variable \x 
e±±£^ rather than v = a 2 , as in hi \m. 

P22P33 ' t~ 1 t-Q l IJ 



Allowing the 23- and 32-entries to be complex conjugates of one another, we further found 
for the corresponding separability function l|, eq. (22)] — where the wide tilde over an i,j 
pair will throughout indicate a complex entry (described by two parameters)— 



n/j, 2 0<fJL<l 

7T /J, > 1 



Further, permitting the 23- and 32-entries to be quaternionic conjugates of one another 



m 



20j, the corresponding separability function eq. (24)] — where the wide hat over an 
i,j pair will throughout indicate a quaternionic entry (described by four parameters) — took 
the form 



5 [5S,M = (Vi S [5S,("»» 2 = (\/i«»('')» 4 ='i; V (3) 

r A* > 1 

So, the real (JTj), complex ((2j), and quaternionic fl3J) HS separability functions accord 
perfectly with the Dyson index sequence (3 = 1,2,4 of random matrix theory {2]]]. "The 
value of j3 is given by the number of independent degrees of freedom per matrix element and 
is determined by the antiunitary symmetries ... It is a concept that originated in Random 
Matrix Theory and is important for the Cartan classification of symmetric spaces 



^ < n < 1 



22 



P- 



480] . The Dyson index corresponds to the multiplicity of ordinary roots, in the terminology 



of symmetric spaces 



23 



Table 2]. However, we remain unaware of any specific line of 
argument using random matrix theory [ijj] that can be used to formally confirm the HS 
separability function Dyson-index-sequence phenomena we have noted above and observed 
in [l|. (The basic difficulty/novelty appears to be that the separability aspect of the problem 
introduces a totally new set of complicated constraints — quartic (biquadratic) in fi jl, eq. 
(5)] eq. (7)] — that the multivariate integration must respect [l|, sec. I.C].) 
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As a further recent exercise, unreported in we found that setting any single one of 
the four components of the quaternionic entry, x 2 3 + W23 + jj^23 + ^23, in the scenario just 
described, to zero, yields the separability function, 



^ < n < 1 



consistent, at least, in terms of the exponent of //, with the Dyson-index pattern previously 
observed. 

Continuing the analysis in [lj], we computed the integrals 

POO 

^sep/ scenario / ^scenario^f^)'^scenarioi,f^)^'f^'i (^) 
JO 

of the products of these separability functions with the corresponding (univariate) marginal 
jacobian functions (which are obtained by integration over diagonal parameters only and 
not any o 
variables 

The ratios of such separable volumes to the HS total volumes 



' the off-diagonal Xy's and yi/s) for the reparameterization of p using the Bloore 

r HS 

sep/ scenario' 



eq. (17)]. This yielded the HS scenario-specific separable volumes V^ s 



00 

yHS = HS I jHS 1 \ r / fi N 

tot/ 'scenario scenario I ^ scenarioKr 1 1 r 1 ^ \ l 

JO 

where cfJ! nario is a scenario-specific constant, gave us in [l| (invariably, it seems, exact) 
separability probabilities. (For the three scenarios listed above, these probabilities were, 
respectively, ||,| and ^.) 

Based on the numerous scenario-specific analyses in [l|, we are led to believe that the 
real, complex and quaternionic separability functions adhere to the Dyson-index pattern for 
general scenarios, when the Hilbert-Schmidt measure has been employed. This apparent ad- 
herence was of central importance in arriving at the assertions in [1, sees. IX. A. 1 and IX. A. 2] 
that the HS separability probabilities of generic [9- dimensional] real and [15-dimensional] 
complex two-qubit states are ^ and Jj, respectively. There we had posited — using mutually 
supporting numerical and theoretical arguments — that jl|, eq. (102)] 

S rea l(p) OC - (3 ~ /i 2 )/i, (7) 

and, further pursuing our basic Dyson-index ansatz (fitting our numerical simulation ex- 



tremely well 



, Fig. 4]), that (S rea i(fi)) 2 oc S comp i ex {[i). (Obviously, we must as well make 



the further claim that (<5yea/(p)) 4 oc S quat (p). Unlike the real and complex cases, however, 
we have performed no numerical analyses such as those in [3] to guide us as to the proper 
coefficient of proportionality to employ. Thus, we have no specific assertion to advance as 
to the two-qubit quaternionic separability probability — although J| or might be readily 
suggested.) 



Now, employing formulas (13) and (14) of Dittmann [24| for the Bures metric — which 
avoid the possibly problematical need for diagonalization of p — we were able to find the Bures 
volume elements for the same three basic (one pair of free off-diagonal entries) scenarios. 
We obtained for the real case, 



-JJ„r,, _ y / 7^TTv / 1 " l>'-\ " p22y / P22 

4a/1 - £23 (P22P 2 + P11) vVp22 + (1 - P11) Pn 



dv Sff = — — ; / = o . ; dp n dp 22 dx 23 dp, (8) 

for the complex case, 



dV2SZ = Pi iP2 2 (pii+P22-l) d Pll dp 22 dx 23 dy 23 dp, (9) 

M Vl-!^-^W + Pn) (-Pi, + Pn + P 2 p| 2 ) 

and for the quaternionic case, 

dV l^% S = gdp 11 dp22dx 2 ady2adu2 3 dv2adii, (10) 



where 



A = -PnP22 (Pll + P22 - 1) 



and 



B = 4 V 1 ~ - V 23 ~ X h - Vis (P22P 2 + Pll) {-p\l + Pll + P 2 P2 2 ) 2 • 

In analyzing the quaternionic case, we transformed — using standard procedures 



m p. 



495] 25!, eq. (17)] — the corresponding 4x4 density matrix into an 8 x 8 density matrix with 



(only) complex entries. To this, we found it most convenient to apply — since its eigenva 



26 1 



27 



ues 
, P- 



and eigenvectors could be explicitly computed — the basic formula of Hiibner 
2664] for the Bures metric. 

Integrating these three volume elements over all the (four, five or seven) variables, while 
enforcing nonnegative definiteness of p, we derived the Bures total (separable and nonsepara- 
ble) volumes for the three scenarios— = £ » 0.822467, V Bur S^ = £ w 0.484473, 

' tot/[{2,3)] 12 ' tot/[( 2 ,3)] 64 

and V Bur S^ = ^ « 0.126835. 

tot/[(2,3)] 768 

We note importantly that the Bures volume elements (@, Qj, (|TU|) ). in these three 
cases, can be factored into products of functions of the off-diagonal variables, u 23 , v 23 , x 23 



and t/23 j and functions of the diagonal variables, pn,p 2 2 and p. Now, we will integrate 
(one may transform to polar and spherical coordinates, as appropriate) just the factors - 

-, 1 = , -, = = and , |== | = — involving the off-diagonal variable(s) over those 
variables. In doing this, we will further enforce (using the recently-incorporated integration- 



over- implicitly-defined-re 



transpose-criterion 



gions feature of Mathematica) the Peres-Horodecki positive-partial- 



101 ]. expressible as 

p 2 -* 2 3 >0 (11) 

in the real case, 

f ~ A - 2/23 > 0, (12) 

in the complex case, and 

^ - ^23 - Viz ~ U \z ~ V 23 > °) ( 13 ) 

in the quaternionic case. (None of the individual diagonal pa's appears explicitly in these 
constraints, due to an attractive property of the Bloore [correlation coefficient /off-diagonal 
scaling] parameterization. Replacing p 2 in these three constraints by simply unity, we obtain 
the non-negative defmiteness constraints on p itself, which we also obviously must enforce.) 
Performing the indicated three integrations, we obtain the Bures separability functions, 



7r p, > 1 

2 sin -1 (p.) < p < 1 



S^(p)={ , (14) 




s mM = \:'r, / (15) 

- p 2 ) < p < 1 



and 



P > 1 



Sp^Qi) = < . (16) 

1 |7r 2 (- v ^?p 2 -2 v ^7?+2) 0<p<l 

Then, utilizing these three separability functions — that is, integrating the products of the 
functions and the corresponding remaining diagonal-variable factors in the Bures volume 
elements ((JS]), (J9])), ( (fit)]) ) over the p,pn and p 22 variables — we obtain separable volumes of 
yJp/[(2, 3 )] = 0.3658435525 and 

V Bures = yBures _ ^2/^ + n L^UQ - 6 + 7r) « 0.124211 (17) 

•ep/[(2,3)] tot/[(2,3)] 32 V ^ 64 V ^ V ^ 




T5 2 » 

3S B JSi e3 {n) 

FIG. 1: Joint plot of the normalized Bures quaternionic separability function — , the square 

$ Bures Qj) 

of the normalized Bures complex separability function [(2 '^ , and the fourth power of the nor- 

malized Bures real separability function [(2, ^ ] . The order of dominance of the three curves is 

the same as the order in which they have been mentioned. 

and consequent separability probabilities, respectively, of 0.4448124200 and (our only exact 



Bures separability probability result in this study (cf. [281])) 



.Bures 4C - 6 + 7T 



pxures = _ w Q.256384, (18) 

aep/[(2,3)] -jr ' V 1 



where C ~ 0.915966 is Catalan's constant (cf. 29]). (This constant appears commonly 
in estimates of combinatorial functions and in certain classes of sums and definite inte- 
grals 30, sec. 1.71.) Further, for the quaternionic case, y BuT fS. — « 0.012954754466, and 

^ ' sep/[(2,3)] ' 

pBure^ w 0.10213883862. (The corresponding HS separability probability was also of the 

sep/[(2,3)] 

same relatively small magnitude, that is, ^ l|, sec. II. A. 3]. We have computed the various 
Bures separable volumes and probabilities to high numerical accuracy, hoping that such 
accuracy may be useful in searches for possible further exact formulas for them.) 

So, the normalized — to equal 1 at [i = 1 — forms of these three separability functions 
are [(2,3)] — ^ [(2,3)i — anc | — vm» — _ i n pjg_ [JJ we plot — motivated by the appearance of 

the Dyson indices in the analyses of [l| — the fourth power of the first (real) of these three 
normalized functions together with the square of the second (complex) function and the 
(untransformed) third (quaternionic) function itself. We find a very close, 

S Bures( \ S B -^ S U) 3S B J^ S (u) 

°[(2,3)] ^ , [(2,3)] K ^\ 2 _ / [(2,3)]^\ / -t q\ 

n ) ~ ( 2vr } ~ ( 4^ ); (19) 

but now not exact fit, as we did find in [lj for their (also normalized) Hilbert- Schmidt 

cffs ( s 5^ (m) 25^ ( M ) 

counterparts M^, JSSfll and (p ? p). 
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As an additional exercise (cf. PJ), we have computed the Bures separability function in 
the case that a single one of the four components of the (2,3)-quaternionic entry is set to 
zero. Then, we have (falling into the same tight cluster in Fig. [TJ when the |-power of its 
normalized form is plotted) 

^ 2 (4-v / 21og(3 + 2 v ^)) p>l 
±tt (fiy/1 - p 2 - sin" 1 ^)) (y/2 log (3 + 2\/2 - 4)) < p < 1 

n 

We have been able, further, using the formulas of Dittmann [24], to compute the Bures 
volume elements for the corresponding (five-dimensional) real and (seven-dimensional) com- 
plex scenarios, in which both the {2,3} and {1,2} entries are allowed to freely vary. But 
these volume elements do not appear, now, to fully factorize into products of functions (as 
is the case for (jSJ) and Qj) involving just pn,P22,P an d just the off-diagonal variables x^-'s 
and yij's. The requisite integrations are, then, more problematical and it seemed impossible 
to obtain a univariate separability function of p. 

For instance, in this regard, we have for the indicated five-dimensional real scenario that 



^fil)!( 2 ,3)] = \\\ BC , D + dp 11 dp22dx 12 dx 23 dfi, (21) 



where 



A = -P 2 U P1 2 (PH + P22 - 1) ((P 2 - 1) P22 + 1) , (22) 

B= (p 22 p 2 + Pll ) 2 ,C = x 2 12 + x 2 23 -l, (23) 

D = (pn + p 22 ) [x\ 2 p 22 {p 22 p 2 + Piif -{{p 2 - 1) P22 + 1) (-Pu + Pn + p 2 p 2 22 )) (24) 
and 

E = -^l 3 p22 (PU + P22 ~ 1) {-Pu + Pll + P 2 pl 2 ) • (25) 

So, no desired factorization is apparent. 

However, the computational situation greatly improves if we let the (1,4) and (2,3)- 
entries be the two free ones. (These entries are the specific ones that are interchanged 
under the operation of partial transposition, so there is a greater evident symmetry in such 
a scenario.) Then, we found that the three Bures volume elements all do factorize into 
products of functions of off-diagonal entries and functions of diagonal entries. We have 

1 



dV\fi u I),(2,3)\ - q\I~~H2 iw 2 TY7 1 TY7 1 rdpudp2 2 dp33dxudx23, 

X {l (Zi4 - 1) 023 - 1) (P22 + P33 - 1) (P22 + P33) 



(26) 



where simply for succinctness, we now show the volume elements before replacing the p 33 
variable by p. (We note that th 
corresponding complex scenario, 



variable by p. (We note that the expression for dV^^L 3 ^ is independent of pn.) For the 



where 



and 



dV H)J^)\ = li % d Pndp22dp 33 dr u dr 23 d6 u de 23 , (27) 



F = -r u r 23 p u p 22 p 33 (pn + P22 + P33 - 1) , (28) 



G = (r 2 u - 1) (r 2 2 3 - 1) (p 22 + p 33 - l) 2 (P22 + p 33 ) 2 , (29) 

and we have now further shifted to polar coordinates, + hjij = r^cos + isin^). 
For the quaternionic scenario, we have (using two sets of hyper spherical coordinates 
{r^efleflef}) and (r 2 3,^g\4?)), 

dV imm = l\j^ d Pndp22dp 33 dr 14 dr 23 d6$d9$ (30) 

where 

F = sin 2 (<) sin (*«>) sin 2 (*g>) sin ) r^p?^ 2 (-pn - p 22 - p 33 + if 2 pf 

(31) 

and 

G = yjl-r 2 14 yjl-r% 3 (p 22 + p 33 - l) 2 (p 22 + p 33 f • (32) 

r Bures 7T 1 



The total Bures volume for the first (real) of these three scenarios is V^mffg. , 2 3 ^ = |j 



0.484473, for the second (complex) scenario, \Z Bm- 55_- = ^- w 0.507339, and for the 

tot/[(l,4),(2,3)] 192 ' 

third (quaternionic), V^" r 52 = « 0.0039119. 

VH h tot/[(l,4),(2,3)] 245760 

In the two corresponding Hilbert-Schmidt (real and complex) analyses we have previously 
reported, we had the results jjj, eq. (28)], 




and [l|, eq. (34)] 



5 [fl 5 4),(2,3)](^) - { ( 33 ) 



, 7r 2 p 2 < a < 1 
W P> 1 
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thus, exhibiting the indicated exact (Dyson sequence) proportionality relation. We now 
found, for the two Bures analogs, that 



<Sfi"4) 6 , S (2,3)] (^) 



7T 



/i = 1 
27T CSC _1 (//) /i > 1 

27rsin" 1 ( / u) < < 1 



(35) 



C Bures ( \ 

[(M),(2,3)]^ 



16tt 2 
16tt 2 



1 



16vr 2 (^1 - 

and, further still, for the quaternionic scenario, 



1 



fi = 1 
/i > 1 

< \i < 1 



(36) 



oBures I \ 
[(M),(^3)] lW 



16tt 4 



8tt 4 (2(^/^31-^)^+^/^!) 



9/i a 



k 9 



7T 



-a/1 



fJ, 2 fX 2 



2a/1 



/! = 1 
yU > 1 
< /U < 1 



(37) 



Employing these several results, we obtained that y Bm ' s 



sep/[(M),(2^)] 



0.096915844, and V Bure J— 

Se p/[(1,4),(2,3)] 



sep/[(i,4),(2,3)] ~ 0.1473885131, 
0.000471134100 giving us real, complex 
and quaternionic separability probabilities of 0.3042243652, 0.19102778 and 0.120436049. 

We see that for values of \i e [0, 1], the normalized forms of these three Bures separability 
functions are identical to the three obtained above ( ffT4"j) . ffl5l) . ffTB"]) ) for the corresponding 
single- nonzero-entry scenarios. While those earlier functions were all constant for fi > 1, we 

1 N 



now have symmetrical behavior about \x = 1 in the form, Sf^^ rio {jj) = Sf^m-io\~ 



In Fig. El we show the analogous plot to Fig. CD, using the normalized (to equal 1 at fj, — 1) 
forms of the three additional Bures separability functions ( (|35|) . (|36|) . (!37|) ). We again, of 
course, observe a very close fit to the type of proportionality relations exactly observed in 
the Hilbert-Schmidt case ((1531. AM]) ). 

We were, further, able to compute the Bures volume element for the i/iree-nonzero-entries 
complex scenario [(1, 2), (1,4), (2, 3)], but it was considerably more complicated in form than 
those reported above, so no additional analytical progress seemed possible. 

Regarding the possible computation of Bures separability functions for the 9-dimensional 
real and 15-dimensional complex two-qubit states, we have found, preliminarily, that the 
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s(aO 




gSBures^ ^ 

FIG. 2: Joint plot of the normalized Bures quaternionic separability function — [(MM 2 . 3)1 — , the 



4 

s B Jsy a (u) 

" 4), (2,3) ' • 
16vr 2 



square of the normalized Bures complex separability function 1(1,4 jg^ 3)1 — , and the fourth power of 
the normalized Bures real separability function [(1,4) ^ ,3)1 — . Over the interval n £ [0, 1], the three 
functions are identical — with the same order of dominance — to those in Fig. [TJ 

corresponding metric tensors (using the Bloore parameterization) decompose into 3x3 and 
6x6, and 3x3 and 12 x 12 blocks, respectively. The 3x3 blocks themselves are identical in 
the two cases, and of precisely the (simple diagonal) form (if we employ hyperspherical coor- 



dinates) that Akhtarshenas found for the Bures metric using the coset par ameterization 



31 



eq. (23)]. They, thus, depend only upon the diagonal entries (while in [31], the dependence, 
quite differently, was upon the eigenvalues). It appears, though, that the determinants — for 
which we presently lack succinct formulas — of the complementary 6x6 and 12 x 12 blocks, 
do depend upon all, diagonal and non-diagonal, parameters. 

The close proximity observed in this study between certain separability results for the 
Hilbert-Schmidt and Bures metrics is perhaps somewhat similar in nature/explanation to 



a form of high similarity also observed in our previous analysis 32]. There, large scale nu- 
merical (quasi-Monte Carlo) analyses strongly suggested that the ratio of Hilbert-Schmidt 
separability probabilities of generic (rank-6) qubit-qutrit states (6x6 density matrices) to 
the separability probabilities of generically minimally degenerate (boundary/rank-5) qubit- 
qutrit states was equal to 2. (This has since been formally confirmed and generalized — in 
terms of positive-partial-transpose-ratios — to arbitrary bipartite systems by Szarek, Bengts- 



son and Zyczkowski in 33[. They found that the set of positive-partial-transpose states 



is "pyramid decomposable" and, hence, is a body of constant height.) Parallel numerical 



ratio estimates also obtained in 



321 ] based on the Bures (and a number of other monotone) 



11 



metrics were also surprisingly close to 2, as well (1.94334 in the Bures case [32l . Table IX]). 
However, no exact value for the Bures qubit-qutrit ratio has ever been established, and our 
separability function results above, might be taken to suggest that the actual Bures ratio 
is not exactly equal to 2, but only quite close to it. (Possibly, in these regards, the Bures 
metric might profitably be considered as some perturbation of the flat Euclidean metric (cf. 



m) 

We plan to continue to study the forms the Bures separability functions take for qubit- 
qubit and qubit-qutrit scenarios, with the hope that we can achieve as much insight into the 
nature of Bures separability probabilities, if not more, than we obtained by examining the 
analogous Hilbert-Schmidt separability functions p|. (In [35[], we had formulated, based on 
extensive numerical evidence, conjectures — involving the silver mean, y/2 — 1 — for the Bures 
[and other monotone metric] separability probabilities of the 15-dimensional convex set of 
[complex] q ubit- q ubi« states, which we wou.d further aspire to test. One may a.so consider 
the use of monotone metrics other than the minimal Bures one [15( — such as the Kubo-Mori 
and Wigner-Yanase.) The analytical challenges to further progress, however, appear quite 
formidable. 
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